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We present a general technique to implement products of 
many qubit operators communicating via a joint harmonic os- 
cillator degree of freedom in a quantum computer. By condi- 
tional displacements and rotations we can implement Hamil- 
tonians which are trigonometric functions of qubit operators. 
With such operators we can effectively implement higher or- 
der gates such as Toffoli gates and C"-NOT gates, and we 
show that the entire Grover search algorithm can be imple- 
mented in a direct way. 

A quantum computer is a device which is capable of 
coherently processing information which is stored in a 
collection of small quantum systems. Much attention has 
been devoted to quantum computers due to the discov- 
ery of algorithms which enable a quantum computer to 
solve certain computational problems much faster than 
any classical computer jl],|^] . In a quantum computer an 
algorithm is represented as a series of unitary operations, 
and with a set of so-called universal gates acting only on 
single two-level systems and on pairs of two-level sys- 
tems, it is possible to produce any unitary evolution on 
the Hilbert space of a collection of two-level systems so 
that any algorithm can be implemented with these gates 
H . This theorem couples the development of quantum 
computing to the theory of classical computing where 
a similar theorem exists, and the complexity of various 
computational tasks has been analyzed simply by count- 
ing the number of universal gates required to perform the 
entire computation. Proposals for practical implemen- 
tation of quantum computing deal with practical issues 
such as identification of quantum systems which can be 
addressed by the experimentalist, but which do not deco- 
here with time, and it is a particularly interesting task to 
find ways to implement the two-bit gates, acting on the 
joint state of a pair of two-level systems (with internal 
states |0) and |1)), or qubits. 

Starting with the ion trap proposal by Cirac and Zollcr 
a number of proposals for quantum computing exist, 
where the individual qubits are coupled to a harmonic 
oscillator degree of freedom, and where two-bit gates are 
implemented by use of the coupling to such a 'data-bus'. 
In the ion trap, the internal electronic or hyperfme states 
of the ions are coupled to the collective vibrational degree 
of freedom due to the recoil during absorption of laser 
light; quantum dots may be localized in an optical cavity 
and communicate via a single mode of the optical field || , 
and it has been proposed to couple Josephson-junction 
qubits by an LC-oscillator mode in an electrical circuit 



In the original ion trap proposal the state of one 
qubit is transferred to the data-bus which is then brought 
into interaction with the second qubit of the gate. In this 
proposal it is essential that the state of the harmonic os- 
cillator is initially cooled to the ground state. In order 
to be able to use an oscillator which is not initially in a 
known state, one can use a scheme which only virtually 
excites it j7|, so that the internal states of the ions are 
completely disentangled from the unknown state of the 
oscillator. It is even possible to use a scheme which dra- 
matically entangles the qubits with the oscillator degrees 
of freedom, to only magically at the end of the opera- 
tion remove all entanglement and produce an effective 
coupling of two qubits which is completely independent 
of the state of the oscillator To produce a unitary 

time evolution of the form exp(i[iAB), where A, B com- 
mute, these proposals use the simple fact that 

exp(iAixA) exp(i\2pB) exp(— iXixA) exp(—i\2pB) 

= exp(-i\ 1 X 2 AB). (1) 

This property can be seen from the Baker-Hausdorf rela- 
tion since A and B commute and the commutator of the 
oscillator position x and momentum p is a constant. Ap- 
plication of (|]J) requires that one can induce interaction 
Hamiltonians proportional to x A and pB. Since these 
operators can be expressed in terms of lowering and rais- 
ing operators such couplings can be induced using the 
resonance condition associated with excitation and de- 
excitation of the oscillator together with the implemen- 
tation of the internal state operators. If A and B act 
on different qubits, we obtain a two-qubit gate. If they 
both involve many qubits, we can produce multi-particle 
entangled states. It has already been shown that if we 
take A = B = J y = J2i \ a yii Eq. (jl]) leads to an effec- 
tive interaction proportional to Jy , which can be used to 
produce a Schrodinger cat like state of the bits P-pT|. 
Throughout this Letter we apply a Pauli spin notation 
for the description of the qubits. The qubit states |0/(1;)) 
are defined as the a z i — —1(1) eigenstates. 

It will be useful to have a geometric picture of the 
contents of Eq. (|l|): Each of the exponential terms on 
the left hand side are displacement operators for the har- 
monic oscillator (conditioned on the internal eigenstates 
of operators A and B) , and the four terms displace the 
system around the rectangular path in Fig. 1 (a). When 
a system is displaced around a closed loop in phase space, 
the state vector acquires a geometric phase factor equal 
to the enclosed area M. In Fig. 1(a) the area is the prod- 
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uct of the lengths of the sides, which due to the operator 
character of these lengths becomes the product of two 
internal state operators, and the resulting phase factor is 
the operator on the right hand side of Eq. (0). 
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FIG. 1. Translations in xp-phase space of the oscilla- 
tor during gate operation: (a) In the evolution described by 
Eq. (^), the oscillator is displaced by the amount X2B along 
the x-axis, then by —XiA along the p-axis, etc., and when 
it eventually ends up in the initial state, a geometric (and 
internal state dependent) phase factor given by the enclosed 
area AiA2^4-B multiplies the state vector of the system, (b) 
By application of an interaction proportional to Cn, the dis- 
placement along the p-direction in part (a) of the figure can be 
rotated into another direction given by the angle 6C, and the 
area enclosed by the solid line becomes A1A2AB cos(9C). To 
perform a Grover search or C"-NOT operation we need to en- 
close several parallelograms with angles which are multipla of 
CO. An effective way to achieve this is to place all subsequent 
parallelograms so that they share one side with the previous 
one, as shown with the dashed curve for the second parallel- 
ogram. With this construction, the multi-bit operation can 
be achieved by traversing only the outline of the combined 
figure. 

The trick contained in Eq. (Q) suffices to produce two- 
bit gates since the operators A and B can be replaced 
by any pair of single particle operators acting on particle 
one and two. The C-NOT gate, which is obtained by 
using A = (cr z i + l)/2, B = a x % and A1A2 = tt/2, can be 
combined with single particle operations to produce any 
unitary operation acting on all the bits || . This method 
in general involves several one and two particle gates to 
produce multi-bit gates. For instance, in Ref. || 4 one- 
bit gates and 3 two-bit gates were used to construct a 
three-bit gate which apart from phase-factors is equiva- 
lent to the C 2 -NOT or Toffoli gate. Experimentally each 
gate corresponds to turning on a given Hamiltonian for a 
certain duration, and therefore each gate adds an exper- 
imental complication and/or possibility of error. In this 
paper we pursue a different strategy for implementing 
multi-bit gates. We will show that one may extend the 
trick in Eq. (|l|) to produce higher order gates directly. 

In we discussed the application of a Hamiltonian 
with continuously varying terms in Ax and Bp, and we 
showed, in particular, that harmonically varying coef- 
ficients on the operators corresponding to bichromatic 
fields, can also be used to produce the operator prod- 
ucts. To extend these result to multi-bit gates we shall 



need a slightly more general interaction which may be 
described by 



H(t) = v(t)Ax + w{t)Bp + r(t)Cn, 



(2) 



where A, B and C are commuting operators acting on 
the internal states of the atoms, n is the number op- 
erator for the harmonic oscillator, and v, w, and r are 
arbitrary functions of time. With this Hamiltonian the 
time dependent Schrodinger equation for the propagator 
idU(t)/dt = H(t)U(t) has the solution 
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(3) 



with 



R(t) = C / r(t')dt 



V(t) = / Av(t') cos{R(t')) - Bw(t') sin(R(t'))dt' 
Jo 

W{t)= I Bw{t')cos(R{t'))+Av(t')sm(R(t'))dt' 



S(t) 



V(t')[Bw(t')cos(R(t'j) 



+ Av(t')sm(R{t'))]dt', (4) 



It is straightforward to check the solution by taking the 
time derivative of U and using the Baker-Hausdorf re- 
lation to simplify the result. In the xp-phase space the 
net action of this propagator is to perform translations 
(x,p) — > (x + W(t),p — V(t)) followed by a rotation by 
an angle R(t) around the origin. Since the functions V, 
W, and R involve the internal state operators A, B, and 
C the translation and rotation is entangled with the in- 
ternal states of the bits. We now generalize the trick 
applied in Eq. (Q) to ensure that V , W, and R vanish 
after a certain time t, such that the harmonic oscillator 
is returned to its initial state, and we are left with an in- 
ternal state evolution operator exp(— iS(r)), where S(t) 
is the area enclosed by the trajectory in the phase space. 
Note, that the expression for S(t) does not involve opera- 
tors referring to the harmonic oscillator, so that the gate 
is insensitive to the initial state of the oscillator. This 
gate can be applied with the oscillator in an unknown 
state, e.g., in a thermal state. 

As a first concrete example of our procedure consider 
three bits which are subject to the time independent 
Hamiltonian 



TT ,U (Jzl+(Jz2 + 1 ( , 1 

H = \l I — x — <j X 3 n + 
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(5) 



where K is an integer. After a duration r = K2n/fl 
the propagator (0) reduces to exp(— in[(a z i + (J z2 + l) 2 — 
l]a x3 /16) = exp(-z7r(er 2l + l)(a z2 + l)a x3 /8), which is 
exactly the Toffoli gate. (We used the fact that for a sin- 
gle particle Pauli operator a, sin(0er) = sin(#)cr). In the 
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ion trap quantum computer the gate can be achieved by 
applying a single pulse of suitably directed and detuned 
fields to the ions. 

The three particle Toffoli-gate, can be constructed so 
easily because the constant term in a z \ + <J Z 2 + 1 can 
be chosen so that this operator squared yields the de- 
sired combination of internal state operators apart from 
a single particle rotation. This technique is not directly 
applicable for more than three particles, and we have 
not been able to devise a similar simple construction 
with only a single Hamiltonian in this situation. In- 
stead we shall produce gates by sequentially applying 
three different Hamiltonians. To make a C" c -NOT op- 
eration, where the first n c qubits control the action of 
the n c + 1 st qubit, we need a projection operator which 
projects into the space where all the n c control bits are 
in the |1) state. Such a projection operator can be ex- 
pressed as a product of single particle projection opera- 
tors ni^iC "*/ + !)A and the C" c -NOT operation may 
be expressed as exp(— in/2 n]=i ^' 2 +1 ^ er xric+ i). The op- 
erators that are easy to make in practice are sums of 
individual particle operators like J z — J = Ym=i ^T"^- 
To turn the sum into a product, we observe that if and 
only if all n c control qubits are in the |1) state, not only 
is the product Yl?=i ^' 2 +1 - > ec l ua l to unity, also the sum 
J z -J = YZ=i vam shes. We now use the Fourier 

transform Y^k=i C0S (2tt^^V) = mS(N mod m) which can 
also be applied to operators so that: 

l k=l v 7 

The C™ c -NOT gate is thus the product of n c + 1 terms 
ex P( 2 (n 7r +i) C0S (^Ti(^ ~ J))vxn c +i) (A = l,2,...,n c + 
!)■ 

To implement a unitary operator which can be writ- 
ten in the form exp(— i/j,Acos(6C)), where A and C 
are internal state operators, we make explicit use of 
the fact that we have an internal state operator ap- 
pearing inside a trigonometric function in the expres- 
sion for 5 (Q). Geometrically, we follow the construc- 
tion of the parallelogram in Fig. 1 (b): First, we ap- 
ply a Hamiltonian proportional to Ap which performs 
a translation along the x-axis. Then a Hamiltonian 
H ~ Cn makes a rotation of the phase space by an an- 
gle 9C: exp(i9Cn)x exp(-iOCn) = cos(8C)x + sin(9C)p, 
and we perform a translation along the p-axis with B 
equal to the identity, etc. The enclosed area is pro- 
portional to Acos(8C) and the propagator has the de- 
sired form exp(— ifiAcos(8C)). By varying the strength 
and duration of the pulses one can control the parame- 
ters 9 and /j,, and using A = o x ,n c +\ and C = J z ~ J 
the parallelogram results in the time evolution operator 
exp(-i/xcos(0(J z - J))<7 x , na+ i). 

By using the operator identity (^J) we can devise a C" c - 
NOT gate by following the outline of n c + 1 such paral- 



lelograms, one after the other. By rotating each parallel- 
ogram, so that the first linear displacement is precisely 
the opposite of the last displacement of the previous par- 
allelogram, we can save half of the translations. Note 
that the sum over I implicit in the J z term in Eq. (j^) 
just amounts to illuminating several qubits instead of a 
single qubit at a time. 

In 1997, Grover presented a search algorithm j|] that 
identifies the single value Xq that fulfills f(xo) = 1 for 
a function f(x) provided, e.g., by an oracle (all other 
arguments lead to vanishing values of the function). If x 
is an integer on the range between and TV — 1 = 2" — 1 , 
the search algorithm is able to find xo after on the order 
of \fN evaluations of the function. Grover's algorithm 
has been demonstrated on NMR few qubit systems [fl2"|| . 
In the following we show how our proposal can be used 
to implement the search algorithm. 

The quantum algorithm first prepares an initial trial 
state vector populating all basis states with equal proba- 
bility. To implement a full Grover search the function 
f(x) has to be a non-trivial function which is imple- 
mented by the quantum computer, but for demonstra- 
tional purposes, the function f(x) can be encoded by 
letting the state of the register undergo a transforma- 
tion where the amplitude of the xq component changes 
sign and all other amplitudes are left unchanged. This 
step can be implemented by writing xq in binary form, 
bobib2...b„~i, and by applying the unitary operator 

^^exp^n( ^ + f (7) 

Below we show how this time evolution may be imple- 
mented with our procedure. Note that the correspond- 
ing effective Hamiltonian vanishes when applied to any 
state where the qubit value (eigenvalue of a z i) does not 
coincide with 26/ — 1, i.e., the state must be the exact 
representation of x$ to acquire the sign change. 

The crucial step in Grover's algorithm is an 'inver- 
sion about the mean', where the amplitude with the sign 
changed will grow in comparison with the other ampli- 
tudes. In the n— qubit computer with N = 2" amplitudes 
c x , the operation can be written c x — > -k ^2 x i c x > — (c x — 
W Sa;' The sum of all amplitudes of the state vec- 
tor can be obtained as any component in the vector 
M\ip), where M is the N x N matrix with unit elements 
in all positions. The inversion about the mean is there- 
fore given by the unitary matrix [Q] 

Ug = - I, (8) 

where / is the N x N identity matrix. 

A straightforward calculation shows that the M ma- 
trix fulfill (sM) k = s k N k ^ 1 AI, and hence we have the 
exponential 

exp ( s M)=J + £^=J+i( e ^-l)M. (9) 
fc=i 
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Thus, by choosing sN = iir, we get exp(sM) = I — 
jfM, which apart from an irrelevant global phase yields 
precisely the inversion about the mean. 

In the standard binary basis, the matrix M couples 
all states to any other state, and it can be written as 
the tensor product H?= (cr x l + 1), where the single qubit 
operators a x i + 1 are 2x2 matrices with unit elements in 
all positions. The inversion about the mean is therefore 
produced directly by the action of the following multi- 
particle operator 



Ug = exp 




(10) 



where we used N = 2". 

Both Uf and Uq can be implemented effectively using 
(^). To implement the function (0), it is easiest to first 
invert all the bits, which have the value zero in xq, so that 
Uf on that state should encode only unit bit values, i.e., 
Uf is precisely the exponential of the projection operator 
in the left hand side of Eq. (^J) . Following the outline of 
the parallelogram in Fig. 1 (b) with A and B equal to the 

we obtain the exponential 



-1 (<r,i 



identity and C = Y^i- 
of one of the terms in the sum on the right hand side, 
and by combining n + 1 such terms one can construct 
the full sum. After application of this simple Uf, the 
qubits encoding the value zero should be flipped back 
again. All qubits should then have their a x components 
rotated into the z-direction, to use again the operation in 
(|J) to implement Uq, which is the same operator, defined 
for the ^-components of the spins. The whole algorithm 
only requires individual access for the single qubit spin 
flips, encoding xq, and for the final readout. An easy 
demonstration experiment where xq = 1111...1 can thus 
be performed without individual access at all (one only 
needs to verify that the number of excited qubits at the 
end equals the total number of qubits). 

In summary, we have presented a technique to pro- 
duce multi-bit gates in quantum computers where all 
qubits are coupled to a joint harmonic oscillator de- 
gree of freedom. We have derived general expressions, 
and we have exemplified the method by an analysis of 
the Grover search and the C ra -NOT gate, which appears 
frequently, e.g., in error correcting codes A re- 

cent preprint JlJ] has addressed the achievements of so- 
called 'concurrent quantum computing', in which access 
to multi-bit interaction Hamiltonians of the form Hia z i 
is assumed. That paper presents ideas for Grover's and 
Shor's algorithm, without suggesting a practical means 
to implement the interaction. Our procedure provides 
a proposal for such implementation: since we can write 

exp(-i^nz cr zi) = exp[-i^cos(7rX)i(o'2; - l)/2)], a sin- 
gle parallelogram as in Fig. 1 (b) suffices to produce this 
operator. 

It is known how to make C 2 -NOT and C 3 -NOT gates 
by means of one- and two-bit gates, but it is difficult to 
make a theoretical comparison of these implementations 



with our proposal, since we build up the desired one-, 
two-, and multi-bit interactions continuously in time. 
From a practical perspective, however, our scheme should 
be really advantageous. The essential operation in the 
Grover search ( |i"o| ) is implemented without access to the 
individual qubits and, e.g., in the ion trap it is much 
easier to implement the Hamiltonian H = ^/( <J zi — l) n 
than just a single term H = (a z i — l)n in that sum. 
In addition, it is an experimental advantage to apply as 
few control Hamiltonians as possible, since imprecision 
in timing accumulates if many operations are needed. 

A feature of our proposal worth emphasizing is that all 
operations are expressed as unitary gates acting on the 
qubit degrees of freedom. The oscillator is certainly im- 
portant, and only at the end of the gates, do the qubits 
actually decouple from the oscillator. One consequence 
is that the initial state of the oscillator does not have to 
be specified. It can be in the ground state, an excited 
state, or even in an incoherent mixture of states, pos- 
sibly entangled with environmental degrees of freedom, 
as long as this entanglement does not evolve during gate 
operation. 
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